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ABSTRACT 

Numerical  solutions  are  presented  for  the  two-dimensional  flow  past 


a  circular  cylinder  in  an  infinite  domain.  The  flow  is  assumed  to  be 
uniform  at  infinity  and  the  cylinder  is  allowed  to  rotate  with  a  constant 


angular  velocity  ft*  * jl  is  chosen  to  be  in  the  range  (0  -  5W/a)  where 
a  is  the  radius  of  the  cylinder  and  W  is  the  mainstream  velocity  at 
infinity.  To  incorporate  visco-elastic  properties  into  the  f low/ an 
implicit  four  constant  Oldroyd  model  is  used,  and  the  resulting  nonlinear 
constitutive  equations  are  solved  in  parallel  with  the  equations  of  motion 
as  a  coupled  set  of  partial  differential  equations.  The  method  of  solu¬ 
tion  used  is  a  finite  difference  technique  with  block  over-relaxation. 

The  results  are  compared  with  those  of  other  numerical  computations  as 
well  as  with  available  experimental  data.  In  particular,  consideration 
is  given  to  the  influence  of  cylinder  rotation  and  of  visco-elasticity 
on  the  drag  and  lift  experienced  by  the  cylinder  and  on  the  streamline 
patterns  and  vorticity  distribution.  _r--- 
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SIGNIFICANCE  AND  EXPLANATION 


For  the  past  three  decades  or  so  rheoloqists  have  attempted  to  derive 
better  and  better  mathematical  models  of  rheoloqically  complex  fluids. 

The  testing  of  such  models,  however,  tor  complex  flow  conditions,  has  to 
some  extent  been  neglected.  This  has  been  partly  due  to  a  lack  of  adequate 
mathematical  techniques  and  sufficiently  powerful  computers. 

This  report  considers  one  model  -  the  Oldroyd  four  constant  model  - 
which  is  able  to  predict  elastic  and  shear-thinning  behaviour  in  a  liquid. 
Numerical  solutions  are  obtained  for  the  two-dimensional  flow  of  such  a 
liquid  past  a  cylinder  which  rotates  about  its  axis.  Such  a  flow  is  of 
interest  because  the  complex  nature  of  the  interaction  of  the  main  flow 
and  the  flow  due  to  the  rotation  of  the  cylinder  represents  a  severe  test 
of  the  ability  of  the  model  to  represent  the  truq  behaviour  of  the  liquid. 


The  resfionsibi lity  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  NUMERICAL.  SIMULATION  OK  NEWTONIAN  AND  V  1 SCO-ELA:  T  U  * 

FLOW  FAST  STATIONARY  AND  ROTATING  CYLINDERS 

Peter  Townsend 

INTRODUCTION 

The  ilenvdt ion  of  more  And  more  sophisticated  rheological  models  of  visco-clast,  i 
Liquids  has  been  the  subject  of  some  considerable  study  over  the  past  three  decades 
or  so,  and  only  somewhat  lesser  attention  has  been  paid  to  using  the  models  in  complex 
flow  situations.  The  reason  for  this  has  to  some  extent  been  absence  of  adequate 
mathematical  techniques,  and  this  is  particularly  the  case  for  one  of  the  early  models 
-  the  Oldroyd  model.  Although  some  analyt ical /numerical  solutions  were  obtained  for 
either  simple  flows  or  with  flow  conditions  severely  restricted  in  some  way  by  simpli- 
flyinq  assumptions,  it  was  not  until  comparatively  recently  11-41  that  progress  was 
made,  anil  solutions  obtained,  for  more  complex  flows  of  an  Oldroyd  model.  Part  of 
this  development  has  been  linked  to  the  considerable  increase  in  power  of  digital 
computers  over  the  last  few  years,  so  that  numerical  models  of  sufficient  accuracy 
are  now  feasible. 

The  main  difference  between  early  semi-ana lytica l/semi-numerica 1  solutions  of 
Oldroyd  models  and  the  more  recent  numerical  solutions  is  that  instead  of  attempting 
to  find  approximate  expressions  for  the  stress  tensor  as  a  function  of  the  velocity 
distributions,  the  stress  components  are  considered  as  dependent  variables  along  with 
the  velocity  components  and  pressure.  In  other  words,  in  addition  to  the  equations  of 
motion  and  continuity,  one  solves  the  constitutive  equations  as  a  set  of  coupled 
partial  differential  equations. 

In  this  paper  we  apply  this  technique  to  a  study  of  two-dimensional  flow  of  an 
Oldroyd  liquid  past  a  circular  cylinder,  the  cylinder  rotating  about  its  axis.  Such 
a  flow  is  of  interest  because  of  its  complexity  and  also  because,  even  for  the  special 
case  of  a  stationary  clinder,  there  is  conflicting  experimental  and  theoretical  data 
in  the  literature  as  to  whether  the  streamline  patterns  are  moved  upstream  or 
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downstream  try  the  presence  ot  elasticity.  The  interaction  of  the  two  basic  flows, 
namely  the  steady  flow  and  the  circulation  due  to  the  rotation  of  the  cylinder,  results 
m  a  strong  positional  dependence  of  the  rate  of  strain  tensor.  The  ability  of  the 
constitutive  equation  to  represent  such  a  complex  flow  is  severely  tested. 

It  has  been  found  in  the  past  I  b]  that  for  some  flows  in  which  two  sepiarate  flow 
mechanisms  interact,  although  the  Oldroyd  model  is  well  able  to  predict  qualitatively 
effects  due  to  shear- th inninq ,  it  fails  to  adequately  represent  elastic  behaviour. 

The  flow  of  a  Newtonian  liquid  past  a  stationary  cylinder  is  certainly  not  a  prob¬ 
lem  which  has  been  neqlected  in  the  literature,  see  for  example  |b  -  111.  For  this 
special  case  we  make  no  attempt  to  add  to  existinq  knowledqe,  but  the  availability  of 
several  sets  of  data  does  provide  a  useful  check  on  the  numerical  method  employed  in 
this  study. 

In  spite  of  some  very  early  experiments  due  to  Prandtl  and  Tietjens  112),  in  which, 
in  addition  to  flow  visualization  studies,  they  also  measured  draq  and  lift  forces  on 
the  cylinder,  there  has  been  comparatively  little  effort  directed  to  the  flow  of  a 
Newtonian  fluid  past  a  rotating  cylinder.  Thoman  and  Szewezyk  113)  have  calculated 
time  dependent  high  speed  flows,  with  Reynolds  numbers  Re  of  30,  200  and  10*’,  and 
Phuoc  Loc  114)  has  presented  results  for  very  much  slower  flows  (Re  <_  20)  and  slow 
rotational  speeds  of  the  cylinder.  Althouqh  the  two  investigations  do  not  overlap, 
they  do  draw  some  contradictory  conclusions.  In  this  paper  we  attempt  to  resolve  some 
of  the  contradictions,  although  since  we  make  the  basic  assumption  of  steady  flow  we 
do  not  attempt  to  consider  the  higher  Reynolds  numbers  considered  by  Thoman  and 
Szewezyk,  and  restrict  attention  to  Re  40.  Experimental  evidence  suggests  that 
above  a  Reynolds  nimber  of  40  the  flow  past  a  stationary  cyclinder  becomes  unstable 
and  one  sees  a  cyclic  time-dependent  flow  pattern  in  which  vortices  are  alternately 
shed  downstream  from  the  top  and  bottom  of  the  cylinder.  Where  possible,  some  com¬ 
parison  is  made  with  the  experimental  data  of  Prandtl  and  Tietjens. 

If  one  introduces  elasticity  into  the  fluid,  then  for  flow  past  a  stationary 
cylinder,  two  basic  questions  have  been  considered  at  lenqth  in  the  literature. 
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Kit  t  1  \  .1v'«  elasticity  havo  .»ny  ot  tort  on  t  hr  r.t  i  rami  me  p.ittetn,  and  it  so  doe*.  it 


shift  t  >  'thMBluw  aj  troam  or  downst  t  ram  h'lativr  ti>  t  h«  cgui  v.»  lent  Mrwtonun 


t  ir.inl  uu's  econdly .  does  the  pr  esence  of  elusti.  ity  int  hmm*  ot  ,K'»n,t  «•  t  he 


Newt  on tan  ^it  .*‘i  on  the  '  S’ 1 1 ride t .  investigation.-  .t  1  '  ask  the  a  me  guest  tons  f »  t 


•  heai  -  thinning  Muvuuu  . 


Prom  t  he  jui  t  e  con*,  idet  able  numhei  of  inwst  nations  carried  out,  see  tot  example 


IIs  -  M  .  it  norms  cleat  that  elasticity  vi«vnMscs  t  hr  Newtonian  dtag  at  low  Reynold* 


tuimhei  •  and  incioase-  the  dtag  at  highet  Reynolds  minbots  aU'vr  unity.  On  tht*  qin*st  u>n 


of  at  team line  shifting,  tu*wever,  thrtr  is  »\vxpiete  disagreement.  I'm*  might  export 


diffeient  f  luid  nvniels  to  give  diffeient  piedictions,  hut  rvrn  expet  imrnt a 1  evidence 


is  cont radiotory .  Perhaps  some*  rmernt  work  duo  to  Mona  124)  may  throw  somr  light  on 


tho  problem  in  that  they  claim  to  bo  able  to  show  experimentally  both  upstteam  and 


downstream  shifting  depending  on  tho  Woissenberg  number.  In  tho  ptesent  wotk,  because 


of  cvMtiputat tonal  teat  net  ions,  wo  are  not  able  to  vary  out  fluid  parametets  suf f it  lent ly 


to  attempt  to  pi  edict  both  upstream  and  downstteam  shifting  although  wo  are  able,  to  a 


ettain  oxtont,  to  isolate  tho  separate  effects  of  elasticity  and  shear-thinning 


behaviour  on  tho  streamline  pattern. 


To  our  knowledge  there  is  no  existing  data,  either  theoretical  or  experimental. 


on  visco-elastic  flow  past  a  rotating  cylindet  although  it  is  hoped  that  some  will  bo 


available  shortly  as  a  result  of  an  experimental  i nvost igat ion  by  R.  H.  bird  and  K.  1. 


Christiansen,  cunently  in  progress  at  the  University  of  Wt  scons  in-Nad t  son .  As  for 


the  stationary  case,  we  again  attempt  to  Isolate  elastic  and  shear -thinning  behaviout , 


and  present  drag  and  lift  coefficients  for  various  Reynolds  numbers,  rotational  speeds 


and  fluid  parameters.  Streamline  pro  lections  and  vorttcity  distributions  are  also 


given  tor  selected  cases. 


ha  n  nu  ky 

We  hsuli>i  thr  t  eady  flow  .»(  a  vii  -•  eU»t  h'  liquid  |  .t*t  a  ii  u!.*i  '  t  f 

«.i  Uu>  4  and  of  mfinlt«  length.  w**  refer  ail  m«t  ion  f<’  .1  wet  .*f  .  \  1  i/utf  J,  .1)  p  la  1 
coot  dmates  U.H.f),  where  the  /  axis  .ohm  1«  ■■  with  t  hr  .ixu  of  the  ylindei  (see 
Figure  1).  At  'flu*  1  a  t  v)  ••  distant  e  away  t  t  *m  the  .  >  1  indet  ,  t  hr  liquid  1  i.  .ij.s'J’H'.I  t.  flow 
left  to  r  iqht  with  a  constant  velocity  W. 

it  we  take  a  velocity  v»»  t.'t  v  11  iu,  v,  n  then  the  relevant  rquat  ions  of  rK^mentum 
and  >ntinuit\  ate 
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wh«-ie  .  1»  the  iensity  of  the  l  i  jui  l,  p  .lenote:>  the  notn'i  u-  pressure  aiul  p 

p  ‘  ,  «nU  y'  the  relevant  >-vun|ionents  of  the  extta  sttess  tensot. 

1  r  “  >  \  t*  “  » 


l  *  r> - 


We  need  now  to  select  a  net  of  v  >n*t  \  t  ut  1  vr  equations  which  will  character  ire  >ui 
v  1  *.  o-e last  1 .  liquid  ind  m  doinvj  wr  would  wish  to  include  both  elastic  behaviour 
and  a  shear  dependent  viscosity.  One  of  the  simpler  models  available  is  the  foui 
constant  Oldroyd  mkIcI  l.'t»]  Although  it  might  Ih»  argued  that  this  model  is  limited 
to  4  iualit.it  ive  picture  of  viscoelastic  behaviour,  it  has  been  used  successfully  in 
the  past  for  a  number  of  different  flow  situation*  nee  l*'4*)).  The  Oldroyd  eguat  ions 
ate  hi  1  wen  by 


'Ik  .  T  'ik 
v  ♦  V  p 

*  1  Tt  1 


*1  auk 

♦  I,  p  e 


„  '  IlHk  .  r  tmk’** 

3n0*  *hpt* 


14) 


wh«r»  n(  i*  the  rero-*heat  viscosity,  V..V,  sn.l  , 

l  •  0 


(U  ik 

are  time  const  a:.ts,  e  1 


n 

the  rate-of-at  t  am  tensor,  and  is  the  co-deformat  tons  l  derivative.  The  relevant 

Pt 

components  of  (4>  are  given  by 


^Parentheses  placed  around  suffi.es  are  used  to  denote  the  physical  oomj'onent  s  ot 

tensors.  The  notation  p  .  .  correspond*  to  -  t  used  in  l^Sl  . 

9  (lk>  1  11 

••Covariant  suffice*  are  written  below,  oont r avar lant  suffice*  above,  and  the  usual 
sunr  %t ion  convention  for  repeated  suffice*  is  implied. 
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-.  -..u.  :»i 


>■«  ►  1 


Tv'  iMini'let*  the  specification  of  the  problem  we  need  appropriate  boundary  condi¬ 
tions  >n  the  ylinder  and  at  infinity.  Since  our  stress  components  are  qiven  by  a 
set  of  differential  equations,  then  in  addition  to  conditions  on  the  velocity  com¬ 
ponent*,  we  requi  t  **  stress  boundary  conditions.  For  the  infinity  boundary  it  is 
omputat lonal ly  convenient  to  specify  a  steady  streaming  in  the  positive  *x’ 
dire  t ion  on  some  very  large  cylinder  t  -  r^,  where  r^  a.  The  conditions  then 
become 
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U  W  v'o  •  ,  V  -  W  I  M' 


*  (I  1  )  *  ( *  tM  (t*H  ) 


n  the  cylinder  w«-  h.ivr  tfi.it 


and  if  the  cylinder  rotate*;.  with  arviul.tr  velocity  then 


where  W  »  a  . 

0 

For  the  stress  component*  w*»  have  no  simple  expulsions  that  we  may  write  town, 
an.t  we  must  look  instead  at  what  form  the  constitutive  equations  take  on  the  cylinder 
surface.  For  brevity  we  restrict  attention  for  illustration  purposes  tv'  the  'rtf' 
component  tvjuat  ion  which  may  be  written  it\  the  form 


•  re  .  fv  j£_^  3u 

p  *  'll?  jT-  -  T7P 


30 


't'C  .>  .  vi 

r  p  ‘  TT  If1 


ftu.v)  m> 


where  f  is  a  Known  function  of  u  and  v.  It  will  be  noted  that  the  only  stress 


omi-onent  derivative  which  appears  in  this  equat ion  is  a  *  t*  •  .tenv.it  ive,  so  that 
provided  the  velocity  distribution  is  known  in  the  tluid.  then  111),  toqether  with  the 
equivalent  'rr‘  and  *  t*  0 '  component  s ,  represent  a  set  of  coupled  ordinary  differ¬ 
ential  equations,  with  cyclic  boundary  conditions,  to  be  solved  fot  the  stiess 
distribution  on  the  cylinder.  This  forms  the  basis  for  an  iterative  scheme  to  be 
described  later. 

before  attempt inq  to  solve  the  system  of  equations  11),  IJ),  IM  ,  tt>)  and  17)  we 
perform  a  ruwbei  of  s  imp  1 1  f  loat  ions .  We  tust  introduci’  a  stream  function  ,  defined 


1  J*  h 

u  *  —  —  ,  v  •  —  — *  . 
r  JO  ’  Jr 

With  this  choice  the  equation  of  continuity  is  satisfied  identically. 

Secondly  we  make  a  1 1 ansform.it  ion  into  non-dimensional  variables  as  follows 


'  T  4  , 


u*  u/W,  V 


v  W 


U  i) 


NUMERICAL  SOLUTION 


When  one  refers  to  the  Oldroyd  model  as  a  'simple*  model  then  one  is  talking  about 
the  'simple'  representation  of  fluid  properties.  The  coupled  nonlinear  partial  differ¬ 
ential  equations  which  one  obtains  for  stress  and  velocity  components  from  using  this 
model  are  certainly  far  from  simple  to  solve,  and  unless  some  very  restrictive  pertur¬ 
bation  analysis  is  employed,  one  must  resort  to  numerical  techniques.  In  an  earlier 
paper  11],  a  numerical  finite  difference  solution  was  presented  for  pulsatile  flow  of 
an  Oldroyd  fluid  in  a  circular  pipe.  In  this  work  the  coupled  stress/velocity  equations 
were  solved  as  an  initial  value  problem,  starting  from  rest,  and  advanced  in  time  until 
a  steady,  or  quasi-steady  state  was  obtained.  In  the  present  paper,  the  problem  is 
treated  as  a  boundary  value  problem  in  which  we  seek  coupled  velocity  and  stress  vari¬ 
ables  consistent  with  conditions  specified  on  closed  boundaries.  This  technique  has 
been  used  by  Walters  and  his  co-workers  f2-4]  for  flow  through  constricted  channels 
and  flow  around  sharp  corners. 

In  order  to  discretize  our  equations  we  need  to  adopt  some  finite  difference  grid. 

A  cylindrical  polar  grid  does  not  give  a  very  satisfactory  distribution  of  grid  points. 
We  need  to  concentrate  grid  points  near  to  the  cylinder  where  variables  are  changing 
more  rapidly,  and  where  greater  accuracy  is  required.  Accordingly  we  set 

5  =  log  r  (20) 

and  work  in  terms  of  (5,0)  coordinates  on  the  logarithmic  grid  shown  in  Figure  2. 
Twenty  grid  points  were  used  in  the  5  direction  and  forty  in  the  0  direction. 
Although  this  is  a  somewhat  coarser  grid  than  perhaps  desirable,  one  is  severely 
limited  by  the  very  considerable  amounts  of  computer  arithmetic  needed  to  solve 
problems  of  this  kind.  5^  was  taken  to  be  3.0  giving  an  outer  non-dimensional  radius 

r*  of  20.1. 

00 

In  terms  of  (5,6)  coordinates  our  modified  momentum  equations  become 

9iji  3m  3i p  3m  1  I- 32m  3^m  . 

[36  35  -  35  36 J  “  «e  [^2 


[  3  /3P(rr) 

3P(r0) 

.  -  -  \  _  _L 

cap 

+  2p 

+  !fkl' 

[30  \  35 

30 

P(rr)  P(66) /  35 

\  35  2P 

r8  30  / 
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and 


With  those  di  t  ference  approx imat ions  the  voiticity  equation  may  !<»•  wiitten  in  t !.« 


E  ^  ♦  E  u>  ♦  E  w  ♦  E  ut  ♦  k  «  *  k 

1  i-lj  2  i  ♦  1  j  3  ij-1  4  i  ]  ♦  1  ('ii  i) 


(.V) 


where  the  K's  And  K,  .  are  functions  of  4/,  P  ( rr )  *  ^  ( rO )  and  0ne  obvious 

solution  technique?  which  can  now  be  applied  to  solve  (26)  is  to  use  successive  ov«*i  - 
relaxation  by  {Hunts  at  each  node  of  the  finite  difference  grid.  In  making  a  sweep 
through  the  grid  the  computer  program  must  incorporate  the  cyclic  boundary  condition 
by  i ecogn i z i ng  t  ha t 

u>.  .  a  u> .  kl,  i  *  1, 2, , . .  ,M  (27) 

1 1  IN 

Attempts  to  apply  such  a  scheme  proved  successful,  but  the  solution  was  extremely  slow 
to  settle  down  to  the  symmetric  form  about  the  x-axis  which  one  expects  in  the  case  of 
a  stationary  cylinder.  A  more  satisfactory  method  which  was  finally  adopted  was  to 
write  (26)  in  the  form 


«,  <k)  .  (k)  Ik) 

E.u)  .  ♦  E_u>  ♦  K-u* .  , 

1  1-lj  5  lj  2  l  ♦  1 3 


B  P  (K)  P  (k-l> 

k  -  E  ,u> .  .  ,  “  E  ,u> .  -  , 

lj  3  13-I  4  ij+l 


(28) 


where  k  is  an  iteration  number.  We  may  then  iterate  in  blocks,  each  block  consisting 
of  a  complete  set  of  u>  values  for  a  particular  '  i'  value.  In  other  words  one 
solves  in  a  series  of  circles  concentric  with  the  cylinder.  The  tri-diagonal  system 
represented  by  (2.8)  is  solved  directly  using  a  modified  Gaussion  elimination 
algorithm  which  incorporates  the  cyclic  boundary  condition.  Essentially  one  solves 
a  system  of  linear  equations  of  the  form 

Au»  *  b  (29) 

where  A  takes  the  form 


A  £ 
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The  discrete  forms  of  the  remaining  stream  function,  and  stress  com|«nent 
equations,  and  also  the  boundary  stress  equations  (c.f.  (Ill),  take  essentially  the  same 
form  as  (28)  and  the  same  block  solution  alqorithm  is  used  in  each  case.  Aqain,  one¬ 
sided  differences  are  used  to  ensure  diaqonal  dominance  for  some  of  the  nonlinear  terms 
in  the  stress  equations. 

The  direction  of  the  block  iteration  used,  varies  accordinq  to  the  boundary 
conditions.  For  the  stream  function  we  iterate  inwards  from  the  non-zero  conditions 
on  the  r^  boundary,  while  for  w  and  p  we  iterate  outwards  away  from  the 

non-zero  conditions  on  the  cylinder. 

In  order  to  obtain  consistent  approximations  for  all  the  variables  an  outer 
iterative  scheme  is  set  up  as  follows 


(0) 


(0) 


(p! 


(0)  (0)  (0)  .  .  . 
rr)P(r8)P<09)  boUndary 


(0)  (0)  (0)11  .  _ 

P(rr)P(r0)P(00)J  lnterl°r 
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(1) 


f  (1) 

[P (rr) 


(1)  (1) 

P  (r@  )P  (68 )  b°Undary 


p!”  p!1’  interior 

r(rr)r  (r@r  (00) 


(31) 


The  outer  iteration  is  continued  until  two  successive  iterates  agree  to  some 
predetermined  tolerance. 

In  order  to  optimize  the  convergence  of  the  scheme  set  out  in  (31),  and  indeed, 
for  most  cases,  to  obtain  convergence  at  all,  some  considerable  attention  must  be  paid 
to  smoothing  both  interior  and  boundary  values  of  variables.  When,  for  example,  a  new 
approximation  i|/n+1'  is  computed,  then  a  weighted  mean  of  this  value  and  the 


previous  iterate  given  by 


7 (n+1) 

* 


,  (n)  ^  ,,  .  , (n+1) 

pji  +  (1  -  p  Jili 
s  s 


0  <  p  <  1 


(32) 


is  actually  used  for  ij/  for  the  next  stage  of  the  outer  iteration  procedure.  For  low 
Reynolds  numbers  and  low  rotational  speeds,  little  or  no  smoothing  was  necessary,  but 
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i.  jji.i  *•»  valu*  iik  r  l  initial  violent  ,  filiation-.  ,r  unbounded  growth  in  t  hi 
v»luti  *n  hal  *  t*»  lumped  out.  K<»t  ; ,  u.  and  the  interior  value;-  of  the  .tresses  a 

•nstant  value  ?  *l  wa  um-.I  for  k  ,  whil<-  for  boundary  stress  and  vorttcity  valu«  •  •  , 

S 

which  |  rove  1  to  t»«  n*>i.  critical,  a  variable  smoothing  wa  .  applied.  Initially  a  i 
value  f  as  high  a  .  ).  *»m  was  necessary  for  some  Reynolds  numbers,  and  this  was  then 

relax*-  1  as  the  iteration  proceeded. 

The  block  iterative  scheme,  in  addition  to  faster  convergence*,  has  apparent  ly 
other  advantages.  Phuo  Loc  (14),  who  used  point  iteration  on  a  similar  grid  to  that 
employed  here  to  solve  the  rotating  cylinder  problem  for  a  Newtonian  liquid,  experienced 
numerical  stability  problems  as  the  rotational  speed  is  increased.  No  such  problems 
were  apparent  for  the  block  iterative  scheme  anil  we  were  able  to  extend  the  l’huoc  Loc 
data  to  a  much  higher  range  of  parameters. 

For  the  vi sco-e last ic  liquid,  however,  convergence  proved  to  be  far  more  difficult 
to  achieve.  If  one  is  prepared  to  juggle  with  smoothing  parameters  then  it  may  well  be 
that  a  solution  is  possible  for  any  physically  reasonable  set  of  fluid  parameters. 
However,  in  our  experience,  the  solution  of  the  full  visco-elast ic  problem  is  very 
expensive  on  computer  time  and  with  a  finite  computer  budget  there  seems  to  be  a  limit 
to  the  degree  of  nonlinearity  in  the  equations  that  one  can  accommodate.  This  limita¬ 
tion  has  been  experienced  by  a  number  of  other  workers  in  this  field  (see  (41).  It 
may  well  be  that  the  type  of  iterative  scheme  outlined  in  this  paper  is  unsuited  for 
hiqhly  nonlinear  problems,  and  that  one  mut»t  turn  to  other  technique's.  The  use  of 
higher  order  collocation  approximations  together  with  more  sophisticated  nonlinear 
equation  solvers  (28)  may  well  prove  to  be  one  possibility. 

The  majority  of  the  computation  of  the  data  presented  in  this  paper  was  carried 
out  on  two  DEC  VAX  computers  with  run  times  varying  from  several  minutes  to  several 
hours  in  extreme  cases. 
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When  a  fluid  flow*  past  the  rotating  cylinder,  the  cylinder  i  ••  subject  t  i  net 
force  m  each  of  the  *x‘  and  ’y"  directions,  namely  the  total  drag  and  total  lilt. 

To  compute  these  we  resolve  the  pressure  and  extra  stress  forces  into  x  and  y 
component s  as  follows 

r  *  -p  CosP  ♦  —  [p*  ,Cost<  -  p*  ..Situ)]  (it) 

x  ^  He  lt(rr)  *(r0)  J 

F  •  -p  Sinti  »  -L  Sint*  ♦  p*  a.  rostl]  ('4' 

y  K  He  lr(rr)  ( r© >  1 

The  stress  components  p*  ,  and  p‘  ..  are  obtained  from  (141  and  the  pressure  distn- 
Mrr)  tr6) 

button  found  by  integrating  the  ’0’  momentum  equation  (2).  If  we  then  integrate  I'M 
and  ('41  around  the  cylinder  surface  we  obtain  non-dimensional  drag  and  lift  coefficient 
and  respectively,  where 


total  drag 


D  P-2a 


total  lift 
pw‘a 


l  IS) 


NPMFKIPAI  KP  PI  V: 


U)  Ncwtonun  t  l  >  *w  j^a •. t  .1  ?.t  at  umui  y  .  }l  mdei 

This  has  been  covered  very  ext  «*ns.  i  ve  1\  in  the  literature  and  w»  do  t,  t 

attempt  to  add  to  what  has  a  1 i  eady  been  said.  Hie  extensive  numerical  4a  t  a  ava  1 1  ab  1  «• 
do  howinwr  offer  one  very  useful  accuracy  check  on  the  work  presented  in  this  paper  . 
particularly  In  view  of  the  somewhat  coarse  finite  difference  grid  used.  In  Figure  » 
w«*  have  plotted  the  drag  coefficient  as  a  function  of  Reynolds  number  .  We  hav* 
rest  r i  c 1  ed  attention  to  a  range  of  Reynolds  numbers  where  one  would  expect  from  experi¬ 
mental  evidence  to  find  a  stable  f low* svwmet t ic  about  the  *x'  axis.  For  higher 
Reynolds  nu ml '♦•is  an  instability  occur*  in  which  vortices  are  periodically  shed  down¬ 
stream  from  alternate  sides  of  the  cylinder,  and  in  this  case  a  time  dependent  analysis 
is  necessary. 

F tom  the  data  given  in  Figure  1  we  see  that  the  present  work  is  in  good  agreement 

with  most  previous  workers  and  we  feel  some  confidence  therefore  in  out  numerical  method. 

It  is  noticeable,  however#  that  one  set  ot  results  due  to  rhoman  and  srewcryk  1101  is 

very  much  out  of  line  with  the  bulk  of  available  data  and  it  is  not  clear  why  this 

should  be.  The  Thorrvin  and  Ssewcsyk  work  does  differ,  to  some  extent,  from  the  other 

invest  igat  ions  in  that  they  do  concentrate  their  work  on  very  high  Reynolds  number  flows 
6 

up  to  Re  •  10  and  to  obtain  numerical  stability  find  it  necessary  tv'  emplos  a  hybrid 
cartesian  polar  grid  system.  What  numerical  effects  this  has  ate  not  clear.  We  only 
point  out  these  differences  between  their  stationary  cylinder  work  and  that  ot  several 
later  invest igat ions  because  of  a  number  ot  differences  between  their  rotating  cylinder 
work  and  the  current  invest igat ion.  Hiese  differences  will  be  appaient  in  the  next 
sect  ion. 

( 1 1 >  Newtonian  flow  pis t  a  rot  a  t l ng  _cj  1 i nde  t 

In  spite  of  some  very  early  flow  vi sual i rat  ion  exper  iments  tor  this  problem  ll-l, 
theorot leal  solutions,  unlike  for  the  stationary  case,  have  been  few  in  n umbel .  The 
only  two  sets  of  numerical  data  of  which  we  are  aware  at  this  time  are  given,  together 
with  the  present  work,  in  Figures  i4'  and  (*>'  where  we  plot  the  drag  and  lift 
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u'cf  t  u  rents,  t  functions  of  tht*  cyl  index xotational  .m  <«d  for  iifferent  Reynold* 
numbers.  It  is  clear  that  Phuts  U>c’s  conclusion  114)  that  the  lr  ag  is  a  slowly 
decreasing  function  of  rotational  speed  ami  that  tht*  lift  is  a  linear  function  of  ret 
tional  speed,  are  somewhat  premature  in  view  or  the  very  limited  i  ange  of  parameters 
considered.  Far  trum  a  decrease,  we  find  a  substantial  increase  in  drag  as  the  cylinder 
is  rotated,  with  some  evidence  that  this  increase  gradually  tails  off  as  tin*  rotat  .  :.al 
speed  increases.  This  result  is  qualitatively  in  agreement  with  the  data  of  Thoman  and 
Szewczyk  (13)  but  some  considerable  difference  is  apparent  in  actual  drag  coefficient 
values. 

From  Figure  5  it  can  be  seen  that  the  present  work  predicts  a  lift  force  which 
rises  initially  to  a  peak  as  the  rotational  speed  is  increased,  but  then  falls  rapidly. 
The  higher  the  Reynolds  numbei  of  the  mainstream  the  sooner  the  lift  experienced  by  the 
cylinders  reaches  a  maximum  and  falls  away.  Again  qualitatively  we  see  similar  behaviour 
in  the  data  of  Thoman  and  Szewczyk  (13)  although  they  predict  very  much  greater  values 
of  the  lift  coefficient.  The  results  of  Phuoe  Loc  (not  shown)  essentially  lie  along 

the  solid  lines  for  Re  **  5  and  20,  but  for  such  a  limited  range  of  K  (W  -  0.5)  that 

0  0  — 

they  offer  no  help  in  resolving  the  discrepancy  between  the  other  two  data  sets. 

In  view  of  the  considerable  change  in  lift  as  the  rotational  speed  is  increased, 
it  is  of  interest  to  plot  the  streamline  projections  for  a  number  of  situations. 

Figures  (6-9)  show  streamline  projections  for  a  Reynolds  number  of  30  and  various 
rotational  speeds.  It  can  be  seen  that  the  pair  of  vortices  formed  behind  the  cylinder 
in  the  stationary  case  are  soon  swept  away  by  the  rotation  of  the. cylinder  and  the 
streamlines  near  the  cylinder  are  drawn  closely  around  it  until  one  sees  a  rinq  of 
trapped  fluid  rotating  with  the  cylinder.  These  results  are  in  excellent  agreement 
with  the  flow  visual izat ion  studies  of  Prandtl  and  Tietjens  (12). 

(in)  Visco-elastic  flow  past  a  stationary  cylinder 

From  the  many  investigations  available  in  the  literature,  both  experimental  and 
theoretical,  two  facts  emerge.  First,  it  seems  generally  agreed  that  the  presence  of 
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elasticity  induces  an  initial  small  decrease  in  drag  at  low  Reynolds  numbers,  which 
charges  to  a  somewhat  larger  increase  in  drag  as  the  Reynolds  number  is  increased. 
Secondly,  there  is  complete  disagreement  as  to  whether  elastic  behaviour  in  the  liquid 
shifts  the  streamline  pattern  upstream  or  downstream  relative  to  the  Newtonian  stream¬ 
lines.  The  various  positions  adopted  by  the  many  investigations  are  set  out  in  Table  1. 
Some  related  work  on  flow  past  a  sphere  is  also  included.  Recent  experimental  data 
presented  by  Mena  (241  suggests  that  either  an  upstream  or  downstream  shift  are  possible 
depending  on  the  value  of  the  Weissenberg  number. 

Of  the  available  theoretical  investigations,  a  number  have  used  implicit  differ¬ 
ential  fluid  models,  but  in  all  cases  very  restrictive  simplifying  approximations  are 
employed  to  make  the  equations  tractible.  In  some  cases,  these  simplifying  assumptions 
have  been  severely  criticised  in  later  work.  The  present  investigation  represents  the 
first  full  numerical  solution  of  the  problem  for  an  implicit  visco-elastic  model  and 
for  the  range  of  parameters  considered,  a  small  downstream  shift  in  the  streamline 
pattern  was  found.  The  results  are  shown  in  Figure  10  where  in  addition  to  the 
Newtonian  streamlines,  those  of  two  other  fluids  are  plotted.  The  equivalent  vorticity 
distributions  are  given  in  Figure  11.  In  an  attempt  to  isolate  the  influence  of  elastic 
and  shear  thinning  properties,  two  sets  of  Oldroyd  parameters  were  selected.  In  one  of 
these  (ig  is  set  to  zero  which  results  in  a  liquid  of  constant  viscosity  but,  provided 
*1  >  ^2’  3  Pos^t^-ve  non-zero  first  normal  stress  difference.  We  shall  make  reference 

in  the  results  to  an  ’elastic'  liquid  meaning  the  constant  viscosity  case,  and  to  a 
'visco-elastic'  liquid  meaning  the  full  four  constant  model  where  shear  thinning  effects 
are  included.  The  results  given  in  Figure  10  indicate  that  the  small  downstream  shift 
due  to  elasticity  is  increased  when  a  shear-rate  dependent  viscosity  is  included.  The 
streamline  projections  for  the  'elastic'  liquid  compare  favourably  with  those  of  Pilate 
and  Crochet  (211  who  used  a  constant  viscosity  second  order  fluid  model.  It  is 
unfortunate  that  computational  limitations  with  the  solution  technique  employed  here 
do  not  allow  us  to  investigate  a  wider  range  of  parameters  and  attempt  to  reproduce 
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TABU:  1:  RESULTS  FHOM  VARIOUS 

AUTHORS.  FOR 

VlSCO-FlAsTlC 

FLOW  PAST  A  STATIONARY  CYLINDER  1 

:m  OR  SPHERE  (») 

Expt 1/ 

Investigators  Theory 

St  i  rami  l  nt's 

Dr  a^ 

Kt 

Leslie  | 1 M  *  T 

Downst.  ream 

Decrease 

1 

Ultman  and  Penn  [lb)'*  T 

Upst ream 

Decrease 

<< 

1 

E 

Downst ream 

Decrease 

<< 

1 

James  and  Acosta  1171*  E 

— 

Increase 

> 

1 

Mena  and  Caswell  1181'*  T 

Downst ream 

Decrease 

« 

i 

Broadbent  and  Mena  119]^*  E 

No  Change 

Decrease 

<  < 

i 

?.ana,  Tie fenbr uck#  E 

and  Leal  120)* 

Upstream 

— 

<< 

i 

Pilate  and  Crochet  )2lS  T 

No  Change 

Decrease 

i 

Downst ream 

Increase 

i 

Siqli  and  Coutanceau  (221*  E 

— 

Decrease 

<  < 

i 

T 

Up st ream 

— 

0 

Ada chi,  Vos hi ok a  E 

and  Sakai  1 211*  TK 

Upst  ream 

No  Shift 

Increase 

Decrease 

(•i 

(•i 

-60) 

-601 

*Non-Newtonian,  inelastic  liquid. 


both  upstream  and  downstream  shifting  suggested  by  the  latest  experimental  evidence 
of  Mena  124). 

We  have  chosen  not  to  present  drag  coefficient  results  at  this  stage  since  they 
will  appear  in  the  next  section  as  a  special  case  of  the  rotating  cylinder  data.  We 
remark  in  passing,  however  that  it  was  found  that  for  small  Reynolds  numbers  (Re  <  0.1), 
elasticity  produced  a  small  decrease  in  drag,  while  for  higher  speed  flows  (Re  >  5)  an 
increased  drag  was  found.  This  observation  again  agrees  with  the  numerical  predictions 
of  Pilate  and  Crochet,  and  with  experimental  data  given  by  James  and  Acosta  [17],  and 
Broadbent  and  Mena  [19).  For  the  parameters  considered  the  shear  thinning  properties 
of  the  liquid  appear  to  have  little  effect  on  the  drag.  It  could  be  argued,  however, 
that  shear-thinning  properties  become  apparent  only  at  higher  shear  rates  or  with 
higher  values  for  anc*  wo  t*1an  considered  here,  and  certainly  we  do  see  later 

a  much  greater  effect  of  shear  thinning  at  higher  rotational  speeds  of  the  cylinder. 

( lv)  Visco-elastic  flow  past  a  rotating  cylinder 

The  flow  past  a  rotating  cylinder  represents  an  interesting  test  problem  for  a 
mathematical  model  of  a  visco-elastic  liquid  in  that  the  interaction  of  the  two 
separate  flow  mechanisms  produces  a  very  complex  flow  situation.  Experience  in  the 
past  [5)  has  indicated  that  for  some  flows  in  which  one  has  interacting  steady  and 
oscillatory  flow  mechanisms,  the  Oldroyd  model  was  qualitatively  well  able  to  predict 
effects  due  to  shear-thinning  behaviour  but  failed  even  qualitatively  to  represent 
elastic  effects. 

Figures  12,  13  and  14  show  drag  and  lift  coefficients  plotted  against  rotational 
speed  for  the  purely  elastic  liquid,  the  full  visco-elastic  liquid  and  the  equivalent 
Newtonian  liquid.  Dealing  with  the  drag  first,  it  can  be  seen  that  elasticity  appears 
to  contribute  an  almost  constant  additional  drag  independent  of  the  rotational  speed. 

The  shear-thinning  behaviour,  on  the  other  hand  has  little  or  no  effect  at  low  rota¬ 
tional  speeds,  but,  as  the  speed  increases,  the  thinning  of  the  liquid  seems  to  reduce 
the  total  drag  below  the  purely  elastic  case,  and  eventually  below  that  of  the 
Newtonian  liquid. 
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The  lift  coefficient  would  to  In*  somewhat  mor  e  sensitive  than  the  dr  ag  to 

the  visco-elastic  properties  of  the  liquid.  At  low  rotational  speeds,  for  both  elastic 
and  visco-e last w  liquid  ,  one  •  .  a  marked  increase  in  lift,  but  eventually  the  lift 

drops  away  to  a  value  below  that  t  the  equivalent  Newtonian  liquid.  The  shear- thinning 
behaviour  of  the  liquid  appear  s  to  accelerate  this  deo.iy  in  the  lift  coefficient. 

In  view  of  the  fairl*  substantial  changes  in  the  forces  on  the  cylinder  brought 
aoout  by  visco-elasticity,  it  is  of  interest  to  see  what  flow  patterns  bring  about 
these  changes.  In  Figures  is  and  lt>  we  plot  streamline  projections  for  two  different 
flow  conditions.  Figures  1'  and  id  give  the  corrvqvr.linq  vorticity  distributions. 

The  first  remark  that  one  an  make  is  that  the  rotat ion  of  the  cylinder  considerably 
amplifies  the  effects  or  vi s.o-elast :icity  competed  with  the  stationary  situation,  and 
one  sees  much  greater  changes  in  the  streamline  patterns.  For  the  lower  rotational 
speed  (Figure  Is*)  one  sees  that  the  elasticity  of  the  liquid  tends  to  draw  the  stream¬ 
lines  in  closer  around  the  cylinder  compared  with  the  equivalent  Newtonian  liquid, 
whereas  the  shear  thinning  properties  have  the  opposite  effect  tending  to  move  the 
streamlines  away  from  the  cylinder.  This  effect  is  exaggerated  as  the  rotational 
speed  increases  (Figure  16),  with  the  result  that  for  the  purely  elastic  liquid  a  much 
larger  ring  of  fluid  rotates  with  the  cylinder  compared  with  the  Newtonian  and  visco¬ 
elastic  cases. 


*This  latter  effect  is  to  be  expected  in  view  of  the  locally  higher  Reynolds  number 
near  the  cylinder  because  of  the  reduced  fluid  viscosity. 
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•  n  thi-  t  <  r  wo  have  derived  a  full  numerical  solution  for  the  flow  of  a  four 
constant  'M  it  yd  liquid  past  a  rotating  cylinder. 

’umparison  of  the  result;  for  a  Newtonian  liquid  and  a  fixed  cylinder,  with  those 
o:  everal  authors  in  the  literature  suggests  that  the  accuracy  of  the  method  employed 

here,  in  spite  of  the  use  of  a  somewhat  coarse  finite  difference  grid,  is  adequate  for 
our  purposes. 

For  the  Newtonian  flow  past  a  rotating  cylinder,  the  rotation  of  the  cylinder  has 
a  significant  effect  on  the  flow,  even  at  low  rotational  speeds,  tending  to  draw  in 
the  streamlines  in  a  'wrap-around'  fashion  around  the  cylinder.  So  much  so  that  at 
higher  rotational  speeds  one  sees  a  body  of  liquid  trapped  in  a  constant  circulation 
around  the  cylinder.  The  drag  and  lift  coefficients  rise  rapidly  with  increasing 
rotational  speed  to  a  peak  value,  and,  in  the  case  of  the  lift  coefficient  one  sees 
an  equally  rapid  decay  as  the  rotational  speed  increases  still  further. 

The  effect  of  elasticity  on  the  streamline  pattern  for  a  stationary  cylinder  is 
to  produce  a  very  small  shift  downstream.  This  shift  is  increased  by  the  presence  of 
shear- thinning  properties.  A  small  decrease  in  drag  is  found  for  both  elastic  and 
visco-elastic  liquids  at  low  Reynolds  numbers,  becoming  a  somewhat  larger  increase  in 
drag  at  higher  Reynolds  numbers. 

When  the  cylinder  is  rotated,  the  elastic  and  shear- thinning  properties  would 
appear  to  have  opposite  effects,  the  former  tending  to  draw  the  streamline  pattern 
in  even  closer  around  the  cylinder  compared  with  the  Newtonian  case,  while  the  shear- 
thinning  properties  seem  to  have  the  reverse  effect.  At  low  rotational  speeds  and 
Reynolds  numbers  in  excess  of  5  an  increased  draq  and  lift  result  from  the  visco¬ 
elastic  properties  of  the  liquid,  but  for  the  shear  thinning  liquid  one  sees  an 
eventual  decay  in  drag  at  higher  rotational  speeds  to  below  the  Newtonian  value.  For 
the  lift  coefficient,  for  both  elastic  and  visco-elastic  liquids,  the  early  increased 
lift  decays  at  higher  rotational  speeds  to  a  value  below  the  Newtonian  lift,  the  decay 
being  somewhat  more  rapid  for  the  visco-elastic  liquid. 
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The  very  expensive  nature  of  the  solution  techniques  used  in  this  invest iqu 


f  »ori 

suggests  that  for  highly  nonlinear  problems  there  is  a  need  for  faster,  more  soph:  ■■■- 
ticated  numerical  algorithms  to  solve  the  coupled  momentum  and  constitutive  equations. 
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